Dynamics of quantum correlations for two-qubit coupled to a spin chain with 

Dzyaloshinskii-Moriya interaction 
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We study the dynamics of quantum discord and entanglement for two spin qubits coupled to a spin chain 
with Dzyaloshinsky-Moriya (DM) interaction. We numerically and analytically investigate the time evolution 
of quantum discord and entanglement for two-qubit initially prepared in a class of X-structure state. In the 
case of evolution from a pure state, quantum correlations decay to zero in a very short time at the critical point 
of the environment. In the case of evolution from a mixed state, It is found that quantum discord may get 
maximized due to the quantum critical behavior of the environment while entanglement vanishes under the 
same condition. Moreover, we observed sudden transition between classical and quantum decoherence when 
single qubit interacts with the environment. The effects of DM interaction on quantum correlations are also 
considered and revealed in the two cases. It can enhance the decay of quantum correlations and its effect on 
quantum correlations can be strengthened by anisotropy parameter. 



I. INTRODUCTION 

Entanglement is thought to be the fundamental resource 
for quantum computation and communication |1]. It is well 
known as a kind of quantum correlation, which is original 
from the superposition principle of quantum mechanics. Re- 
cently, quantum discord is realized as a different kind of quan- 
tum correlation other than entanglement, it is arising from the 
difference between quantum mutual information and maxi- 
mum of quantum conditional entropy [2, 3]. This nonentan- 
glement coiTelation may be applied to the quantum compu- 
tation, which is based on deterministic quantum computation 
with one quantum bit (DQCl) protocol 10. HI. Moreover, A. 
Brodutch et al. |6] have shown that the changes in quantum 
discord is pointed to be an indicator of failure of a local oper- 
ations and classical communications (LOCC) implementation 
of the quantum gates. Both entanglement and quantum dis- 
cord characterize the same quantumness features of quantum 
system in the pure state, namely, quantum discord reduces to 
be entropy of entanglement. However, they become two dif- 
ferent measures for quantum system in the mixed state. For 
instance^quantum discord may be nonzero in the disentangled 
state \^M- 

The properties of quantum discord are intensively investi- 
gated due to its potential application. It is realized that quan- 
tum discord may be generated in the Heisenberg models at 
finite temperature [9-13]. Quantum discord behaves in a very 
different way from that of entanglement by changing the pa- 
rameters. For instance, the situation where quantum discord 
increases with temperature while entanglement decreases is 
revealed lfl3ll . Entanglement is also known as an indicator 
to signal a quantum phase transition (QPT) fl4[|, so does the 
quantum discord. Sarandy has analyzed the quantum discord 
in the critical systems, such as XXZ and transverse field Ising 
models, and shown that quantum discord exhibit signature of 
the QPT ifTsll . Furthermore, Werlang et al. have show that 
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quantum discord, in contrast to entanglement, spotlight the 
critical points associated with QPT for Heisengber model even 
at finite temperature 1 16]. More recently. It is found that quan- 
tum discord detect quantum critical points associated with 
first- and higher-order QPTs caused by field and three-spin 
interactions at finite temperature ITtIi . 

Decoherence of the quantum system due to interacting with 
its environment is the main obstacle to practice quantum com- 
putation tasks. Thus, it is necessary to reveal the dynami- 
cal properties of quantum discord for designing the protocol 
to against decoherence. The dynamics of quantum discord 
both in the Markovian and nonMarkovian environment are 
theoretically and experimentally investigated fls'-'Jl']. The 
researchers have found the situation where quantum discord 
disappears in the asymptotic time limit but entanglement un- 
dergoes a sudden death lis'], which is found in the Morko- 
vian environment. In this sense, quantum discord is thought 
to be more robust than entanglement. Recently, there is an 
increasing investigation on decoherence due to spin environ- 
ment |22-28], such as single qubit coupled to the environment 
and two qubits coupled to the environment. It is revealed that 
quantum coherence is dramatically destroyed in the critical 
point of the QPT from the external environment. 

In this paper, we investigate the dynamical behavior of 
quantum discord and entanglement for two-qubit coupled to a 
spin chain with DM interaction. The DM interaction is arising 
from the spin-orbit coupling |29], which often appears in the 
models of low-dimensional magnetic materials. The dynam- 
ics of the system also depends on the initial state. We choose 
the two-qubit to be a class of X-structure state, which is gen- 
eral to contain Bell-diagonal state and Werner state, which 
are very important and famous in the quantum information 
theory. We start in section by introducing and diagonaliz- 
ing the model, and we calculate the reduced density matrix 
of two-qubit. In section |III] we analytically and numerically 
evaluate quantum discord and entanglement, and present the 
main results and discussion. The last section is devoted to the 
conclusions. 
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II. MODEL AND ITS SOLUTION 

We consider the environment to be an general XY chain 
with z-component DM interaction. The total Hamiltonian for 
two spin qubits transversely coupled to the chain is given by 



with 
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H = + H, 
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(1) 



•C^/Cr/+l)' 



A+6 , 1-6 



(2) 
(3) 



where H^^^ is the self-Hamiltonian of the environment, Hj de- 
scribes the interaction between two-qubit and environment. 
cr^ia = x,y,z) is the Pauli matrix on l-th sites of the chain. 
(T^^gj describes the two-qubit. y measures anisotropy of ex- 
change interaction in the XY plane, A is the strength of trans- 
verse magnetic field, D is the strength of z-component DM in- 
teraction. g-^(g^) describes the coupling strength between 
qubit A{B) and surrounding chain. The parameter 6 controls 
anisotropy of coupling strength of qubit with its environment. 

is the total sites of XY chain. The periodic boundary condi- 
tions cr^_|_i = cr^ is assumed. The eigenstates of the operator 



(^cr^ + ^CT^) are simply given by 

I01> = |00>,|<^2> = |O1),|03> = |10>,|<^4> = 111), 



(4) 



where |0) and |1) denote spin up and down respectively. In 
terms of these states, the Hamiltonian ^ can be rewritten as 



H 



/i=i 



where the parameters Af, are 



(5) 



(6) 



To 



/li(4) = A + g, A2(3) = A + g6, 

and //^'''' is obtained from Z/^''* by replacing A with /l^. 
investigate the dynamics of two-qubit, we need to compute 
the time evolution operator U{t) = exp(-iHt). For the pur- 
pose, we map the projected Hamiltonian //^''' into a one- 
dimensional spinless fermion system with creation and anni- 
hilation operators cj and c/ via Jordan- Wigner transformation 



0-, = 



Y\il-2clc,,)(c, + ci), 

S<1 

-iY](l-2clc,)(c,-c]), 



s<l 

1 - 2c,' c,. 



(7) 



Following by Fourier transforms of the fermionic operator 
given by 4 = 2f ae-'^""'!'^ with k = -M, . . . , M and 
M - (N - l)/2, the Hamiltonian is transformed to momentum 
space, and then using Bogoliubov transformation 



m,A^ = COS -^dk - / sin -^d\. 



(8) 



with e) 
Hamiltonian 



(•y sin — \ 
— — , finally we get the diagonalized 
cos ^ J 



(9) 



where the spectrum A^'*''' is A*''''' 



2(e, 



2£) sin 2^) with 



= ^(i^-cosM)2 + ^2sin2M. 



Suppose the initial state of total system is disentangled with 
PtotiQ) = Pab(O) ® P£(0), where pab{0) and pe{0) are the ini- 
tial state of two-qubit system and environment, respectively. 
P£(0) = |i/'£(0))(^^£(0)| is assumed to be the groundstate of the 
environment. The evolution of the total system is governed by 
Profit) - U{t)p,ot{Q)U\t). Then the reduced density matrix of 
two-qubit AB is obtained by tracing out the environment 

PABit) = TrE[ptot{t)] 
4 

= YjPl^-^tWpABm4>v)\cl>^;}{<py\ (10) 



With F^At) = {'^E\Uf'\t)uf\t)\ilfE), Where uf\t) = 
expC-///^"'?) is the projected time evolution operator driven 

Now we suppose that the two-qubit are initially prepared in 
the X-structure states, 

PAb(O) = \i^lAB + Z ® (^i) (1 1) 



with Iab is the identity operator on two-qubit system, a - 
X, y, z, and the parameters c„ are choose to be real parameters 
that promise the pab(O) be a legal quantum state. This state 
is general to contain Bell-diagonal states and Werner states. 
According to Eq. ( fTOl l and (fTTT l. the reduced density matrix 
can be written in the standard basis {|00), |01), 1 10), 1 1 1)) 



(12) 



with F = {cx - Cy)F\A,{f) and Q — (c^ + Cy)F23(t). To evalu- 
ate the quantum discord in this density matrix, it is necessary 
to compute the decoherence factor F^v(0- Let \G},\ and \G}a^, 
denote the groundstates of the self-Hamiltonian //^^ and the 
projected-Hamiltonian, respectively. By using the transforma- 
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FIG. 1: Quantum discord (solid line) and entanglement(dash-dot 
line) as a function of time t for three different A. Other parameters 
are set as 7 = 1, g = 0.05, S = 0,D = 0, and N = 600. (a) A = 0.85, 
(b)i= 1.00, (c)A= 1.15. 
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FIG. 2: Quantum discord (solid line) and entanglement (dash-dot 
line) as a function of time t under different D. Other parameters are 
set as 7 = l,A = l,g = 0.05,(5 = 0, and A' = 600, (a) D = 0.0, (b) 
Z) = 0.1,(c)Z) = 0.5, (d)Z)= 1.0. 



side measurement 



with 0('*''> = (Ol^"^ - ef)/2 and followinga tedious calcula- 
tion, we obtain the decoherence factor as 11251] 

M 

\F,y(t)\ - Yl^l- sm\2e[''^) sin^CA*-*"'?) - sin^isf'^) 

k>0 

X sin2(A*'*''^f) + 2 sin(20^''"^) sin(20f •'^) 

X sin(A^'^''^f) sm(A[^'h) cos(A^'^''^f - A^ "'f) - 4 

X sin(20^''''^) sin(20^'*''^) sin2(0*'*"' - 

M 

X sin2(A<'''''f) sin2(A<'''>f)]'^' = ^(0- (13) 

k>Q 

When \F^y(t)\ — > 1, it reveals that the two-qubit is slightly 
disturbed by environment. When \Fijy{t)\ — > 0, the two qubtis 
system undergoes strong decoherence due to the environment. 
We present the detail results and discussion in the next section. 



III. CORRELATIONS FOR THE TWO-QUBIT 

Quantum discord is proposed to measure quantumness cor- 
relation in the bipartite quantum system, quantum discord for 
two-qubit state pab is defined as 

Q(Pab) ^ I(Pab) - C(pab), (14) 

where I(pab) = S(pa) + S{pb) - S{pab) with S(p) = 
-Tr{plog2p) is the quantum mutual information, which mea- 
sures the total correlation in state pab- Pa{B) is the reduced 
density matrix of pab- C(pab) in the above definition is the 
maximum of quantum conditional entropy by performing one 



CipAB) = max{5(pA) - J] P/5(p®)}, (15) 

where {11^) denotes a set of von Neumann measurement on 
B, p^* = T rBi^'gPAB^'^)! Pi is the state of A after obtaining 
outcome / on B, where p, - T rABiJ^'gPAB^'g)- C(pab) is sug- 
gested to measure classical correlation in the state pab- 

To obtain the quantum discord of pABit), we need to calcu- 
late the quantum mutual information and classical correlation. 
The quantum mutual information is directly calculated, 

4 

I[PAB(t)] = 2 -I- ^ W„ log2 (Dn, (16) 

n=l 

with wi = \(\ - c, + |Q|), W2 = 5(1 - c- - (jj-i = 
\(\ -H Q -i- |r|), and W4 - \(\ -HQ - |r|). Notice that the condi- 
tion S [pa(0] - S [pBit)] satisfied, we have the unique value of 
classical correlation irrespective of whether performing mea- 
surement on subsystem A or B. In order to compute the 
quantum discord for Eq (fT2l l. we propose the complete set of 
von Neumann measurement of subsystem B, {Tl^g\i - 1,2) = 
{|<A<i))<0(i>|, |<*(2))<;^(2)|)^ with 10(1)) ^ cos6l|0> + e'^sineil) and 
10(2)) ^ e-'Vsin6/|0) - cos 011), the parameters 6,ip e [0,27r]. 
Hence, we have the remain state of subsystem A with an out- 
come /, 

(0 _ f 5[1 - (-l)'c: cos(26l)] ^-^(e-'fa + e'-^T) sin(20) 
Pa - \^ (-irl(^^i<p^* ^ g-'Vp) sin(20) i[l -H (-l)'c, cos(26l)] 

(17) 

and the respective probability p, = j. Consequently, we ob- 
tain the classical correlation of Eq. ( [12] ) 

ClpABit)] = ^ log2(l ^ log2(l - §), (18) 
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where & = maxjlc^l, '"^'^' j. Finally, the quantum discord is 
given by 

Q[pAB(t)] = ^[(l-C, + |Q|)log2(l-C, + |Q|) ; 

+(l-c,-|Q|)log2(l-c,-|Q|) 

+(i + c,-|r|)iog2(i+c,-|r|) 

+(1 + c, + in) iog2(i + c, + in)] - cipABit)] 

(19) 

Entanglement is well known as a kind of quantum corre- 
lation, which is different from quantum discord. In order to 
investigate the relation between quantum discord and entan- 
glement, we choose the entanglement of formation (E) for 
quantifying the amount of entanglement of two-qubit. The » 
£ is a monotonically increasing function of concurrence (C) »» 

iH, 

E = -f(C') log2 f(C') - [1 - f(C')] log2[l - f(C')], (20) 

with f(C') = (1 + Vl - C'2)/2. The concurrence C is defined 
as C = max{0,2^,„„ - Tr[^jpABPAB]], where pab = ct^ ® 
cTgpABtr^ ® cTg, ^max is the maximum eigenvalue of ^JpabPab- 

One can directly calculate concuiTence for the state given by FIG. 3: (a) Quantum discord and (b) entanglement as a function of 
( , |„i , 1 time t and X. Other parameters are set as y = 1, g = 0.05, 6 = 0,D = 

Eq. (inil, C = max |0, ^-^^j- 0, and N = 600. 



A. evolution from pure state 



To investigated time evolution of quantum correlations of 
the two-qubit initially prepared in the pure state, we set the 
parameters d = 1, and -c,. - Cz - 1, then the initial state 
becomes Bell state -^(|00) -i- |11)). Substituting these values 
into Eq. ( fT9b , and noticing that )? = 1, i.e., classical correla- 
tion always equals to 1, finally, we obtain the quantum discord 
for this case, Q = HMill iog2(l - \Fu{f)\) + h]MlA \og^(i + 
|Fi4(f)|). Obviously, Q is a monotonically increasing function 
of variable |Fi4(f)|, gets its maximum (2=1) with |Fi4(f)| = 1 
and minimum {Q - 0) with |Fi4(f)| = 0. In a previous work. 
Yuan et al. revealed \F{t)\ decay to zero in short time evolution 
when A approaches the critical point Ac - 1 |26] in the weak 
coupling regime {g <K 1). Thus one can expect that Q just 
vanishes at the critical point. To evaluate the entanglement, 
we need to calculate the concurrence and it is C - |/^i4(f)l 
for this case. From above discussion, both quantum discord 
and entanglement are the monotonically increasing function 
of variable |Fi4(f)|. Thus, entanglement behaves in a similar 
way as quantum discord does. In Fig.[Tl we plot the time evo- 
lution of quantum discord and entanglement for Ising chain 
(r = 1 ), and other parameters ^ = 0.05 , 5 = 0, D = 0, = 600 
at different points A - 0.85, 1.00, 1.15. Quantum discord is al- 
ways less than entanglement in the process of evolution. At 
point A - \, quantum discord and entanglement jointly decay 
to zero at the same time. 

Now, we consider the effects of DM interaction on quan- 
tum correlation. First, let us check the dynamical property of 



correlation at critical point. One may recall the approxima- 
tion of |Fi4(f)l given in Ref.fz?]. Following the similar pro- 
cedure, one can introduce a cutoff number Kc and define the 
partial product for |Fi4(f)l, |i^i4(f)lc = Uto^kit) > \Fu{t)\. 
The partial product can be written as Sit) - ln|Fc(OI = 

-I^lli ln|^^t(OI- One has A^-*' * |^ - 1| + 4Dsin^^j and 
K^^"^ ~\A^,-\\ + 4sin(^j, (ju = 1,4) for small k and large 
A^. Thus one has sin(20^''''') = p^^^^- (j" = 1,4) and 
sin(0*'''* - O*-*-**) * |(^,I^xif-i)l ■ o™i"ing the terms 



related to the sum of k'^IN'^, one obtain the approximation 
for the partial product S (t) 



1 4;ry^g- 

2 A'-(^-l): 



,2( sin-(Al-'^'> 



TJiV + ^^KtmHr '^os(4gO[. Fmally, when A 1, 



e-(T,+T,)/^ with Ti 
= 2567:^y^gD/iA 



in the short time one has |Fi4(f)| ~ 
327rW/('i - lfi:f'k^/N^ and 
1)^ Zf' k^/N^- Consequently, decay of quantum discord and 
entanglement may be enhanced by introducing DM interac- 
tion. The numerical calculation plotted in Fig. |2] shows the 
effects of DM interaction, which is consistent with theoreti- 
cal computing. It is found that the smaller values of y, the 
stronger the effect of DM interaction on the decay of the deco- 
herence factor |27[1 . It implies that the effect of DM interaction 
on quantum correlation may be controlled by the parameter y. 
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FIG. 4: Quantum discord (solid line) and entanglement (dash-dot 
line) as a function of time t and different S. Other parameters are set 
as 7 = l,g = 0.05, = 0, and A' = 600, (a) 6 = 0.0, (b) 5 = 0.1, (c) 
.5 = 0.5, (d) 5= 1.0. 



FIG. 5: Quantum discord (solid line), classical correlation (dot line), 
and total correlation (dash line) as a function of time t. Other param- 
eters are set as 7 = 1,/1 = l,g = 0.05,(5 = 1.0,D = 0, and A' = 600. 



B. evolution from mixed state 

In this subsection, we study the time evolution of corre- 
lations when two-qubit initially prepared in the mixed state. 
We choose the parameters Cx - I, -c,, = c, = 0.2, thus the 
initial state becomes a mixed state Pab(O) - 0.6\^^}{^^\ + 
0.4|^+)<>I'+|, where |<1)+) = ^(I00> + |11)) and = 

-^(|01) + |10», and it is easy to check that Eq. (fTTI) is a fam- 
ily of Bell-diagonal state for Cj = 1, c,. = -Cz,c- e [0,1]. 
To illustrate the dynamical properties of quantum discord and 
entanglement, we carry out the numerical calculation from the 
exact expression (fT9T l and ( l20l i. In Fig. [3] quantum discord and 
entanglement are plotted as a function of time t and A with the 
parameters g = 0.05,5 = 0,i=l,D = 0,A^ = 600. It is ob- 
served that quantum discord gets maximized at critical point 
but entanglement rapidly decays to zero and vanishes in the 
area where quantum discord has nonzero values. This phe- 
nomenon is also different from the former case of pure state. 
It implies environment may enhance quantum discord of two- 
qubit. 

Let us turn to study the effects of the parameters 6 on the 
quantum correlation. For the case of 6 - 0, the decoherence 
factor 1/^23(01 = 1, otherwise one has 1^23(01 ^ 1, and partic- 
ularly, 1^23(01 = \Fi4(t)\ when 5 = 1. In Fig.|4l;a)-(d), we plot 
the quantum discord and entanglement against time t with 6 
equals to 0.0,0.1,0.5, 1 .0, respectively. It is observed that the 
decay of quantum discord can be enhanced by increasing S, 
however, entanglement seems to be no sensitive to the change 
of 6. When 6 - I, it means one qubit interacts with the sur- 
rounding chain and the other is free from the environment. 
Recalling that |Fi4(f)| decreases to zero in a very short time 
under some reasonable condition, and so does \F23{t)\. Thus, 
the nondiagonal elements in expression ( fT2] i vanish all at once, 
that is to say the coherence of two-qubit quickly disappears. 
As plotted in the Fig.|4jd), quantum discord remains constant 
for an interval of time, and then decays to zero. This situation 
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FIG. 6: (a) and (c) plots of quantum discord as a function of time t 
under different D, and (b) and (d) plots of entanglement as a function 
of time / under different D. Other parameters are set to g = 0.05, 6 = 
0.5, andA' = 600, (a) and (b) 7 = 0.8, (c) and (d) 7 = 0.2. 



seems to be sudden transition between classical and quantum 
decoherence, which has been detail discussed in the Ref. ll32ll . 
Figure. |5] shows that the situation is indeed the transition. In 
the left side of dash-dot line in the plot, one can observe that 
classical correlation decreases and quantum discord maintain 
constant. Whereas, in the other side of the line, classical cor- 
relation do not change with time and quantum discord start 
decreasing. 

To end this subsection, let us focus on how the DM interac- 
tion affects the quantum correlation. In Fig. |6] we plot the 
time evolution of quantum discord and entanglement when 
two-qubit coupled to a general XY chain under different D. It 
is observed that the decay of quantum discord can be slightly 
enhanced by decreasing values of y for case of D = 0. How- 
ever, we find that the decay of entanglement can be slightly 
increased by increasing y, which is opposed from that of quan- 
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turn discord. Nevertheless, entanglement quickly vanishes in 
a much shorter time than quantum discord does. Fig. |6] also 
shows that the effects of DM interaction on quantum discord 
is more remarkable than its on entanglement. 

IV. CONCLUSIONS 

In summary, we have studied the time evolution of quan- 
tum discord and entanglement for two-qubit coupled to a spin 
chain with DM term. We evaluated the quantum discord 
and entanglement for two-qubit to be prepared in a class of 
X-structure state. We have separately studied the two-qubit 
evolutes from pure state and mixed state. The difference of 
quantum discord and entanglement becomes drastic in the 
case of mixed state. In the case of pure state, it is found that 
quantum correlation rapidly decays to zero at critical point 
where environment has a QPT. Moreover, the DM interaction 
can enhance the decay of quantum correlation for the case. In- 
terestingly, quantum discord may get maximized in the case of 
mixed state when environment at the critical point while en- 



tanglement vanishes when environment approaches the criti- 
cal point. We also considered that qubit coupled to the en- 
vironment with different coupling strength from each other's, 
which is controlled by the parameter 6. It has been shown that 
the coherence of two-qubit rapidly vanishes when single qubit 
interacts with the environment. Besides, we have observed 
sudden transition between classical and quantum decoherence 
for the case. The effect of DM interaction in this case is not 
only enhancing the decay of quantum discord and entangle- 
ment, but also enhancing the increasing of quantum discord. 
Furthermore, the effect of DM interaction on quantum discord 
can be strengthened by the anisotropy parameter. However, 
the DM interaction has slight effect on entanglement with re- 
spect to that of quantum discord. 
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